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Abstract: The evolution of net-proton rapidity spectra with ^s^n hi heavy relativistic 
systems is proposed as an indicator for local equilibration and longitudinal expansion. 
In a Relativistic Diffusion Model, bell-shaped distributions in central collisions at AGS 
energies and double-humped nonequilibrium spectra at SPS show pronounced longitudinal 
collective expansion when compared to the available data. The broad midrapidity valley 
recently discovered at RHIC in central Au + Au collisions at ^snn= 200 GeV indicates 
rapid local equilibration which is most likely due to deconfmement, and fast longitudinal 
expansion of the locally equilibrated subsystem. A prediction is made for Au + Au at 
v /5^= 62.4 GeV. 
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Net-baryon rapidity distributions have proven to be sensitive indicators for local equi- 
libration and deconfinement in relativistic heavy- ion collisions PjJ |2] . Whereas statistical 
model analyses of multiplicity ratios of produced particles |H| appear to be consistent with 
the assumption that the system reaches equilibrium and can therefore be described by a 
temperature and a chemical potential, it is clear that one has to consider the distribution 
functions of the relevant observables in order to determine whether the system has indeed 
reached, or gone through, thermal equilibrium. 

Distribution functions of produced particles in transverse momentum or energy have 
been shown in many analyses to be close to thermal equilibrium, if one takes into account 
the collective transverse expansion of the system jl]. Here, the expansion velocities rise 
as the energy increases from AGS via SPS to RHIC when they are estimated based on 
thermal distribution functions. 

In longitudinal direction, however, it is more difficult to determine the degree of equi- 
libration, and the collective expansion velocity. It is the main purpose of this note to 
provide a schematic macroscopic approach that allows to obtain both quantities. There 
have been earlier attempts to extract longitudinal collective velocities at AGS and SPS 
energies by assuming that the rapidity distributions are thermal, and that the remaining 
large difference to the data is then due to expansion [5]. However, the net-baryon rapid- 
ity distributions at AGS, SPS and RHIC energies are clearly nonequilibrium distributions 
even in central collisions |2j, and also the longitudinal distributions of produced particles 
are not fully thermalized |16j . 

To account for the nonequilibrium behavior of the system, the evolution of net-proton 
rapidity spectra with incident energy ^snn= 4.9 to 200 GeV is studied analytically in a 
Relativistic Diffusion Model (RDM) 0, and compared to AGS SPS [Zj and RHIC 
data. In addition to the nonequilibrium-statistical evolution as described in the RDM, 
collective longitudinal expansion is considered. 

The model El El E] is based on a generalized Fokker-Planck equation (FPE) for 
the distribution function R(y,t) in rapidity space, y = 0.5 • ln((E +p)/(E — p)), 



with the nonextensivity-parameter q ^2], and the rapidity diffusion coefficient D y that 
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contains the microscopic physics, and accounts for the broadening of the distribution 
functions through interactions and particle creations. 

The FPE can be solved analytically in the linear model case q=l, with constant D y , 
a linear drift function 

J(y) = (Veq-y)/T v , (2) 

and the rapidity relaxation time r y . This is the so-called Uhlenbeck-Ornstein process 
for the relativistic invariant rapidity 111 j - The equilibrium value is y eq = in the 
center-of-mass for symmetric systems, whereas y eq is calculated from the given masses 
and momenta for asymmetric systems. 

Using 5— function initial conditions at the beam rapidities for net baryons, analytical 
solutions of Eq.((TJ) are obtained for various values of T int /Ty (Fig.l bottom as an example 
for Pb + Pb at SPS). Here the interaction time Ti nt is the time from nuclear contact to 
freeze-out, or the integration time of Eq.Q. It determines how close to equilibrium the 
system can come, and it is obtained from dynamical models or from parametrizations 
of two-particle correlation measurements. In central Au + Au collisions at 200 A GeV, 
this yields Ti nt ~ 10 fm/c jTHj- For known interaction times, the rapidity relaxation times 
and diffusion coefficients can then be obtained from the data. Otherwise, T int /T y and the 
rapidity width coefficient 

r„= [8 ■ ln(2) ■ D y ■ Ty} 1 / 2 (3) 

can be determined from the peak positions and width of the data, respectively, through 
the use of the equations for the mean values and the variances [2] that are derived from 
Eq.flU). 

The transport coefficients D y and r y are, however, macroscopically related to each 
other through a dissipation-fluctuation theorem with the equilibrium temperature T. This 
relation has been derived in from the requirement that the stationary solution of 
Eq.(£Q) is equated with a Gaussian approximation to the thermal equilibrium distribution 
in y-space. At fixed incident energy, this weak-coupling rapidity diffusion coefficient turns 
out to be proportional ^1] to the equilibrium temperature T as in the analysis of Brownian 
motion (Einstein relation) 

D y oc -. (4) 
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When compared to heavy-ion collision data at energies above the SIS region of 1-2 GeV, 
rapidity spectra that are calculated using this relation are consistently narrower Jl] than 
the data because collective processes (in particular, collective longitudinal expansion) 
are not included in the weak-coupling dissipation-fluctuation theorem. Hence, the width 
coefficients T y as obtained from Eqs. ()3l4|) are replaced by effective values T^f. Depending 
on the system and the incident energy, they are typically factors of 2-5 larger than the 
calculated values. 

In this work, I assume that the whole discrepancy Y e Jf — T y is due to longitudinal 
expansion. I obtain the mean collective expansion velocity in indirection for baryons of 
rest mass m and relativistic mass m from the relativistic velocity 



In this expression, the mean energy content of the expansion is written as E co u = //2 • 
(T e ff — T) with an effective temperature T e ff and the equilibrium temperature T in 
analogy to the classical expression for the mean energy. Hence, the expansion enhances 
the particle mass m to its relativistic value m. 

The thermal equilibrium temperature T should be associated with kinetic freeze-out. 
Typical values of T at RHIC energies are 110 MeV. This is significantly below the chemical 
freeze-out temperature of about 170 MeV [B| that is obtained from fits of hadron abun- 
dancies. In the weak-coupling case without collective effects, T e ff=T such that v co u = 0, 
and the distribution functions remain too narrow as compared to the data. With expan- 
sion, the diffusion coefficient is enlarged to its effective value, and since it is proportional 
to the temperature, T e ff is obtained from T using the same enhancement factor. 

With Eqs. (j5l6J) . the collective velocity becomes 




(5) 



For f degrees of freedom (f=l in rapidity space) the effective mass is 



m = m + -■ (T eff - T). 



(6) 
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Comparing the solutions of Eq.Q to Au + Au central collision (5 per cent of the cross 
section) data at AGS-energies y / s^N= 4.9 GeV 6j, it turns out that due to the relatively 
long interaction time r^j and hence, the large ratio Ti nt /T y ~ 1.08 (Table I), the system is 
in rapidity space very close to thermal equilibrium, with longitudinal collective expansion 
at 1^=0.49 (Fig.l, upper frame). The bell-shaped experimental distribution is in good 
agreement with the solution of Eq^ The distributions remain bell-shaped also at lower 
energies jO]. 

This situation changes at the higher SPS energy of ■ s /snn= 17.3 GeV. Here net-proton 
Pb + Pb rapidity spectra corrected for hyperon feeddown [7] show two pronounced peaks 
in central collisions, which arise from the penetration of the incident baryons through the 
system. The gradual slow-down and broadening is described using Eq. ^ as a hadronic 
diffusion process with subsequent collective expansion, v^ oll =0.75 (Table I). The associated 
nonequilibrium solutions with expansion for various values of T int /T y are shown in the lower 
frame of Fig.l. The system clearly does not reach the dash-dotted equilibrium solution. 
Hence, both nonequilibrium properties, and collective expansion are required to interprete 
the broad rapidity spectra seen at the SPS. 

Within the current framework, no indication for deconfmement of the incident baryons 
or other unusual processes can be deduced from the net-proton rapidity data at AGS and 
SPS energies, because the crucial midrapidity region is here too small. 

This is, however, different at RHIC energies ^snn= 200 GeV. The RDM nonequilib- 
rium solution exhibits pronounced penetration peaks with collective longitudinal expan- 
sion 1^=0.93, but it fails to reproduce the BRAHMS net-proton data 8 in the broad 
midrapidity valley (solid curve in Fig.2, bottom): The diffusion of the incident baryons 
due to soft scatterings is not strong enough to explain the net baryon density in the central 
rapidity region. The individual nonequilibrium solutions R\ and i?2 are Gaussians, and if 
they fit the data points near y=2, they necessarily grossly underpredict the midrapidity 
yield because it is in their tails. 

This central region can only be reached if a fraction of the system undergoes a fast 
transition to local thermal equilibrium, dashed curve in Fig.2, bottom. With collec- 
tive expansion of this locally equilibrated subsystem of 22 net protons (t>j,' o ^=0.93), the 
flat midrapidity BRAHMS data are well reproduced in an incoherent superposition of 



5 



nonequilibrium and equilibrium solutions of EqQ 



dN(y, t = Ti nt ) 
dy 



NtRtiy, r int ) + N 2 R 2 (y, r mt ) + N eq R l e ° g c (y). 
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The fast transition of a subsystem of N eq ~55 baryons, or 22 protons to local thermal equi- 
librium suggests that the associated participant partons have reached local equilibrium 
through the sudden enhancement in the number of degrees of freedom that accompanies 
deconfinement. Microscopically, the large gap to midrapidity is thus bridged through 
hard scatterings of partons with subsequent thermal equilibration, rather than diffusion 
of nucleons. The fact that even hard partons can participate significantly in equilibration 
processes is evidenced by the high-pr suppression found in Au + Au at RHIC. 

The local equilibrium distribution R l ° q c with expansion is a solution of Eq.Q for time 
to infinity with an enlarged (effective) diffusion coefficient that is related to the effective 
equilibrium temperature. The nonequilibrium solutions Ri,R 2 of Eq.fJTJ) are calculated 
with the same value of the diffusion coefficient. The underlying picture is that expansion 
affects nonequilibrium and equilibrium distributions in a similar fashion. 

In a schematic calculation for the lower RHIC energy of a/stvv= 62.4 GeV, I have 
used the rapidity diffusion coefficient from 200 GeV to obtain the result in the upper 
frame of Fig.2. Two pronounced penetration peaks can be seen, together with a narrow 
midrapidity valley. The corresponding data have been taken, and are presently being 
analyzed by the BRAHMS collaboration Once they are available, an adjustment of 
and Ng° c may be required. One may expect another correction because the total net 
proton number could change, although the baryon number is conserved. Note, however, 
that this change was found to be negligeable at SPS energies. 

The RDM-results can be further refined if the linear drift function is replaced by 



with the transverse mass m± = ym 2 + because this yields exactly the Boltzmann 
distribution as the stationary solution of Eq^for q=l. This entails, however, to solve Eq^ 
numerically. A numerical solution is also required if one tries to account for multiparticle 
effects through nonextensive statistics |T2| with an explicit nonlinearity in the FPE, 1 < 
q < 1.5. In this case, an approximate result may be obtained from a linear superposition 



a ■ m±sinh(y) 



(9) 
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of power-law solutions of Eq^ EH QJ. The local equilibrium fraction is 

<."°c[l- (1 - q)^cosh(y)}^ (10) 

+oo 

with J i^'dy = 1. For ^s^7= 200 GeV, m±= 1.1 GeV, T = 110 MeV, q=1.4 and 

— oo 

N e? =22 protons, EqJTU] yields the dotted curve in Fig. 2 (bottom). It lies in between 
the Boltzmann-Gibbs local equilibrium results without and with expansion. Hence, it 
simulates to some extent the collective effects. 

Whereas local equilibration in rapidity space occurs only for a small fraction of the 
participant baryons in central collisions at RHIC energies, the opposite is true for produced 
hadrons. This can be inferred from recent applications of the RDM with three sources 
(located at the beam rapidities, and at midrapidity) and ^-function initial conditions [16 . 
About 78% of the produced hadrons in central Au + Au collisions at 200 A GeV are found 
to be locally equilibrated in pseudorapidity space, whereas the corresponding number is 
17% in central d + Au. 

To conclude, I have interpreted recent results for central collisions of heavy systems at 
AGS, SPS and RHIC energies in a Relativistic Diffusion Model (RDM) for multiparticle 
interactions based on the interplay of nonequilibrium and local equilibrium ("thermal") 
solutions. In the linear version of the model, analytical results for the rapidity distri- 
bution of net protons in central collisions have been obtained and compared to data. 
The enhancement of the diffusion in rapidity space as opposed to the expectation from 
the weak-coupling dissipation-fluctuation theorem has been interpreted as collective ex- 
pansion, and longitudinal expansion velocities have been determined from a comparison 
between RDM-results and data based on a relativistic expression for the collective velocity. 
A prediction for net-proton rapidity distributions at y / Siv r /v= 62.4 GeV yields a smaller 
midrapidity valley than at 200 GeV which will soon be compared with forthcoming data. 
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TABLE I. Parameters for heavy relativistic systems at AGS, SPS and RHIC-energies. 
The beam rapidity is expressed in the cm. system. The fit parameter Ti nt /T y determines 
how fast the net-baryon system equilibrates in rapidity space. The second free parameter 
is the effective rapidity width coefficient Ty* . It includes the effect of expansion. The 
deduced longitudinal expansion velocity is v^ oll . At 62.4 GeV, T e Jf will need adjustment(*) 
to forthcoming data. 



Lab System ^/sTn Vb" 1 ' r int /r y T e y ff v^ oU 
(GeV) 

AGS Au + Au 4.9 1.60 1.08 1.45 0.49 

SPS Pb + Pb 17.3 2.91 0.81 2.43 0.75 

RH1C Au + Au 62.4 4.20 0.34 2.94* 0.86 

RH1C Au + Au 200 5.36 0.26 4.29 0.93 
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Figure captions 



FIG. 1: Net-proton rapidity spectra in the Relativistic Diffusion Model (RDM), solid 
curves. Central Au + Au at AGS energies v /s/v r /v = 4.9 GeV (data from [6J is close 
to a longitudinally expanding equilibrium distribution (top). Central Pb + Pb at 
SPS energies ^snn = 17.3 GeV (data from j2j) remains relatively far from thermal 
equilibrium and therefore, double- humped with two penetration peaks (middle). 
Dashed curves are thermal equilibrium distributions without collective expansion, 
dash-dotted curves with expansion. The bottom panel shows analytical Pb + Pb 
RDM-solutions for r mt /r y = 0.08, 0.1, 0.2, 0.5, 0.8, 1.0, 1.1, 1.3, all with longitudinal 
collective expansion, Table I. The transition from bell-shaped (AGS) to double- 
humped (SPS) is clearly shown in the RDM. 

FIG. 2: Net-proton rapidity spectra in the RDM (solid curves) for central collisions of 
Au + Au at RHIC energies y/s^jy = 62.4 GeV (top), and 200 GeV (middle, central 
5% of the cross section) compared to BRAHMS data 8J at the higher energy. The 
net proton content is 158. The high midrapidity yield is not attainable through 
soft scatterings in the Relativistic Diffusion Model. Instead, this region containing 
~ 22 protons at 200 GeV is populated through hard scatterings with ensuing local 
equilibration (dashed curve), and collective expansion (dashed areas). See Table I 
for parameters. The dotted curve is the nonlinear local equilibrium solution of the 
FPE for q=1.4, cf. text. 
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